Revisiting the flocking transition using active spins 
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We consider an active Ising model in which spins both diffuse and align on lattice in one and two dimen- 
sions. The diffusion is biased so that plus or minus spins hop preferably to the left or to the right, which 
generates a flocking transition at low temperature/high density. We construct a coarse-grained description 
of the model that predicts this transition to be a first-order liquid-gas transition in the temperature-density 
ensemble, with a critical density sent to infinity. In this first-order phase transition, the magnetization is 
proportional to the liquid fraction and thus varies continuously. Using microscopic simulations, we show 
that this theoretical prediction holds in 2d whereas the fluctuations alter the transition in Id, preventing 
for instance any spontaneous symmetry breaking. 



PACS numbers: 87.18.Gh, 05.65.+b, 45.70.Vn 

Active matter systems are driven out-of -equilibrium by 
the injection of energy at the single particle level ||T|-|4). 
This microscopic breakdown of detailed-balance is respon- 
sible for a wide range of phenomena that have aroused the 
interest of physicists, from bacterial ratchets HHD to self- 
propelled clusters ll9 tiTTj to the appearance of long-range 
order in 2d systems with continuous symmetries lfT2l[T3l . 
Furthermore, much of this rich phenomenology is captured 
by simple models. For instance, the patterns found in high 
density motility assays could be accounted for using sim- 
ple flocking models lfl4l [T31 while clustering in bacterial 
mixture was successfully modelled using self-propelled 
rods ifTol . 

Nevertheless, despite the successful description of many 
experimental phenomena, a clear-cut understanding of the 
underlying mechanisms sometimes remain elusive. For in- 
stance, even though the flocking transition is one of the 
central features of active matter, it remains one of the most 
debated questions in the field. In their seminal work, Vic- 
sek and co-workers [12 | showed that self-propelled parti- 
cles that align locally can exhibit a transition to long-range 
order in 2d. Initially thought to be continuous lfl2l . this 
transition was later shown to be first order using large scale 
simulations and a finite-size-scaling akin to that of mag- 
netic phase transitions ifTTl . Many works were also devoted 
to nematic lfT8W20ll or metric-free interactions ET1 . the lat- 
ter yielding a continuous transition ll22l . Related flocking 
models were also studied in Id ||23ll24l . where the transi- 
tion was found to be continuous, casting even more confu- 
sion in the field. Beyond the presence of strong finite-size 
effects, a major difficulty in obtaining conclusive numer- 
ical evidence comes from the lack of a theoretical frame- 
work to analyse the finite-size scalings of flocking models. 

In parallel to the numerical studies, much effort was thus 
devoted to the construction of such an analytical descrip- 
tion of the flocking transition. Phenomenological continu- 
ous equations were originally used to prove the existence 
of long-range order in 2d lfl3t . Constraining these equa- 
tions to correctly capture the physics of microscopic mod- 
els is however difficult and an increasing number of studies 



aimed at deriving coarse-grained description starting from 
microscopic models l25] - |29l . While the Vicsek model is 
probably the simplest to simulate, it is one of the hardest 
to coarse-grain, being defined off -lattice, in discrete time 
and involving many-body interactions. Many approaches 
thus relied on simplified models fl25l even though progress 
was recently made on the original problem [26 1. Lots 
of effort was again devoted to the nematic case E7H291 
or the topologic interactions l28l |3T1 . In most cases, the 
coarse-grained equations are quite complex and their ana- 
lytical solutions difficult. Most studies of these equations 
were thus restricted to the linear stability analysis of homo- 
geneous solutions or the comparison between simulations 
of coarse-grained equations and microscopic models. Re- 
cently, non-linear profiles could be computed analytically 
for self-propelled rods [28 1 but all in all, despite the im- 
portant progresses made during the last few years, a unify- 
ing theoretical framework of the flocking transition is still 
missing. 

We present below a tentative step in this direction 
through the introduction of a microscopic lattice model 
which turns out to be much simpler both to simulate and 
describe analytically than its off-lattice counterparts. This 
allows us to bridge micro and macro, and show that the 
flocking transition can here be seen as a traditional liquid- 
gas transition in the temperature-density ensemble, with a 
critical density sent to infinity. 

Let us consider a one-dimensional lattice of L sites on 
which particles have Ising spins s = ±1. Each particle 
hops asymmetrically at rate D(l + se) and D(l — se) to 
its right and left neighboring site. (In higher dimensions, 
the hoping rates are chosen symmetric in all directions but 
one.) There is no exclusion between particles and we note 
nf the numbers of ± spins on site i so that the local den- 
sities and magnetizations are given by pi = + and 
rrii = nf — r\ . The particles also align their spins: on 
site i a spin s changes sign with rate exp(— s(3 — ) where 
/3 = 1 /T plays the role of an inverse temperature. When 
D = 0, the system thus amounts to L d independent mean- 
field Ising models. When D > and e / 0, three dif- 
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FIG. 1. Top-left: Phase diagram in 2d. The red and blue lines 
correspond to low and high densities of phase separated profiles; 
they enclose the region where such profiles can be seen. D = 1, 
e = 0.9, L = 300. Bottom: Snapshots of the different profiles 
average over the transverse direction. Top-right: Phase diagram 
predicted by the RMFM. Red and blue lines correspond to p h and 
pi while black and green dashed lines correspond to the spinodal 
lines where homogeneous profile lose linear stability. D — v = 
r = 1. 



ferent configurations are typically found in the system (see 
Fig. [TJ: at low temperature a uniform ordered phase with 
a non-zero magnetization, at high temperature a uniform 
disordered phase, and in between a phase- separated system 
with high density ordered bands (p, ~ p h , rrii ~ m h ^ 0) 
connected through narrow interfaces to a disordered homo- 
geneous background (pi ~ p e , rrii ~ 0). The stability of 
these profiles in the thermodynamic limit depends on the 
number of spatial dimensions but they are all long-lived in 
finite systems. Let us now show how a simple theoretical 
framework can be constructed to account for the phase di- 
agram of Fig. [T] 

Many coarse-graining approaches used in the past rely 
on factorization approximation of microscopic kinetic 
equations fl2"5ll27ll35ll36ll . On a Id lattice, this amounts to 
a simple mean-field approximation: f((nf}) = (/(n 4 )). 
While this approximation may be quantitatively wrong, it 
often captures phase diagrams of lattice-gas models ex- 
actly [33 1 even in quite complicated cases 11341 . Intro- 
ducing continuous variables x = i/L, v = 2De/L and 
D = D/L 2 and taking the large L approximation, one 
finds for the coarse-grained average fields p{x) = (p^ and 
m(x) = (rrii) the following mean-field equations 



p = Dd xx p — vd x m 



f3m 



m = Dd xx m — vd x p + 2p sinh 2m cosh 

P P 



(1) 

(3m 



(2) 



In higher dimensions, one simply replaces d xx by a Lapla- 
cian A and we use this more general form from now on. 



Looking for the onset of a flocking transition, we lin- 
earize the dynamics for m <C p, which yields 



m 



DAm - vd x p + 2m(j3 - 1) 



m 



a- 



(3) 



where a = f3 2 (l - (3/3) [37]. The solution p{x) = 
p , m(x) = is thus linearly unstable for f3 > 1. Further- 
more, simulations of Eqs. (fTJ and §2§ show that for f3 > 1 
clusters are never stable [38]: they always spread to reach 
the homogeneous ordered profile m(x) = m . Mean- 
field thus predicts a continuous transition from m = to 
m = m (f3) at f3 c = 1, in clear contradictions with Fig.[T| 
This approximation is indeed only valid for p —> oo; for 
finite densities, we need to include 1/p corrections to the 
parameters of the Eq. (T][3 1 and in particular to the critical 
temperature j3 c = 1 + ^, leading to a modified Eq.|3j 

m = DAm — vd x p + 2m(j3 — 1 ) — a— — (4) 

P P 2 

The phase diagram corresponding to Eqs. ([T]) and Q, 
which form our refined mean-field model (RMFM), is pre- 
sented in the top-right corner of Fig. [T] When T < 1, 
homogeneous disordered (resp. ordered) profiles are al- 
ways linearly stable at low enough density p < p ± (resp. 
high enough density p Q > p 2 ). Since p ± < p 2 , there is 
a finite intermediate region p 2 ] where neither homoge- 
neous profiles are stable. In this region, the system sep- 
arates in two homogeneous phases connected with sharp 
fronts: a disordered region with low density pg < pi and 
an ordered region with high density p h > p 2 and m h ^ 0. 

Propagating shocks can be computed analytically when 
f3 is close to 1 by linearizing Eq. Q around the density 
pi = r/{(3 — 1) at which the homogeneous disordered 
profile becomes linearly unstable. We first solve Eq. ([TJ, by 
neglecting the diffusion term in a reference frame moving 
at speed c, to get p as a function of m: 



p(r) = p e + -m(r) 

c 

Eqs. Q and Q then yields for m 
DAm+c(l 



v 2 .^ . v TO 3 

— )d x m+p[p l -p 1 -\ — m(x)\m-a—2 
c 1 c pi 



(5) 







(6) 
0) and 



one gets 



where p = 2r/p\. Looking for ascending (q + > 
descending front (q~ < 0) solutions 

m(r) = ~2^[1 + ta.nh(q ± x)] (7) 

mh\foL 4r 4r 

- m h = — Pe = Pi - — 

9a 

(8) 

Such solutions are consistent with our approximations 
since < 1 and DAp < vd x p when j3 ->• 1. In 
this regime, simulations of the RMFM and Eqs. (T5]{8]) yield 
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FIG. 2. Left: Fraction of the ordered liquid phase when den- 
sity po is either increased or decreased between 1.4 < pi and 
4 > ph for the RFMF (top) and in microscopic simulations (bot- 
tom). Right Corresponding profiles of the system. Parameters: 
RMFM L = 100, v = D = 1, r = 1.6, /3 = 1.75, Ap = 10~ 2 
every At = 15 000; 2d lattice model L = 250, /3 = 2, D = 1, 
e = 0.9, Ap = 10~ 2 every At = 500. 

the same profiles and band velocities. For larger /3, the 
DAp term makes front and rear interfaces asymmetric and 
c > v : the flocks fly faster than the birds. 

Since p e , p h and m h do not depend on p , increasing the 
average density at fixed temperature simply increases the 
fraction of the high-density region. In the thermodynamic 
limit, phase separated profiles can be seen from p e to ph- 
One always has p e < pi < p 2 < Ph so that the clusters 
and the homogeneous profiles are both linearly stable in the 
intervals [p t ,pi] and [p 2 , p h ]- 

The refined mean-field scenario thus resembles an equi- 
librium liquid-gas phase transition in the temperature- 
density ensemble, the total magnetization being propor- 
tional to the fraction of the liquid phase. Varying the den- 
sity po at fixed T, one indeed observes the traditional hys- 
teresis loops shown in Fig.|2j Increasing p , homogeneous 
disordered profiles are seen up to pi where a discontinu- 
ous jump takes the system into a phase-separated profile. 
A further density increase results in a widening of the liq- 
uid phase which almost fill the system when p ;$ p h . (The 
finite widths of the interfaces connecting p £ and p h prevent 
phase-separated profiles for p close to pi/h in finite sys- 
tems.) Going down, the homogeneous liquid phase remains 
metastable until p = p 2 and discontinuously jumps to a 
coexistent state. The fraction of gas then increases until it 
fills the system at p ~ pi . 

Unlike equilibrium liquid-gas transitions, dense and di- 
lute phases in flocking models have different symmetries 
due to the coupling between density and orientation. One 
thus cannot circumvent the transition and continuously 
transform the system from a gas to a liquid: the transition 
line cannot stop at a finite point in the (T, po) plane and, 
indeed, the critical density is infinite J39j. 

Let us now turn to simulations of the 2d active Ising 
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FIG. 3. Reversal of a Id cluster due to a localized fluctuation. 
v 2 is greater than v\ until p(x) = p^ in the whole cluster. (See 
movies in 133.) p = 5, D = 1, e = 0.9, /9 = 1.7. 



model. Beyond the structure of the phase diagram (see 
Fig. [TJ, the RMFM correctly captures the mechanism of 
the transition. The coexistence between homogeneous and 
phase-separated profiles is confirmed and changing p at 
fixed f3 inside the coexistence region simply changes the 
fraction of the liquid phase (see Fig.[2]and l42ll "). 

The scenario presented here can be related to the mea- 
surement of the binder cumulant G = 1 — ,,J done 
in the literature ifTTl [29). Since the high density bands 
have a minimal size £ c , the apparition of a flock in a 
finite-size system corresponds to a non-zero magnetiza- 
tion m ~ m h £ c /L. When this happen P{m) thus has 
a three-peak structures around m = (supercooled gas 
phase) and m = ±m (phase-separated profiles) whose 
relative weights vary across the transition. This leads to 
a negative peak in G which is thus the signature of the 
coexistence between phase- separated profiles and homoge- 
neous disordered profiles. (Another peak could in princi- 
ple be seen close to p h for the coexistence with ordered 
homogeneous profile.) Assuming a sum of three Gaus- 
sians of variance a, the value of the minimum of G is 
G m = -[12(a/m ) 2 + 36(<x/m ) 4 ]~\ which is only 
markedly negative when m S> a. Since both m and a 
vanish when L — > oo, the finite-size scaling of the binder 
cumulant in the temperature-density ensemble is difficult: 
whether the negative peak of G grows IfTTl or recedes 11291 
may simply depend on whether a or m vanishes faster 
and receding negative peaks thus do not necessarily mean 
a continuous transition. Last, many studies use noise rather 
than density as a control parameter, which makes the finite- 
size scaling even more difficult since p e , p h , m h — and 
whence m — also depends on T. 

Let us now consider the one-dimensional case where 
fluctuations strongly alter the transition. First, all three 
profiles (homogeneous ordered, disordered, and phase- 
separated) exist and are linearly stable in finite sys- 
tems 00). The general scenario predicted by the refined 
model thus holds: homogeneous profiles between pi(T) 
and p 2 (T) are linearly unstable and tend to phase- separate 
between linearly stable low-density disordered and high- 
density ordered regions. 

The impact of fluctuations on the transition can be un- 
derstood by looking at the flipping of a whole domain. In 
Id, this happens frequently since an excess of, say, positive 
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FIG. 4. Left & Center: P(jmj) for p„ = 4; Xm{p)/L; (3 = 
1.538, D = 1, e = 0.9. Right: Cluster length as a function 
of time, showing a linear spreading between reversals. D — 1 

e = 0.9 P = 2 po = 3. 



spins on a single site suffice to flip an approaching negative 
cluster (see Fig. [3]). The 2d counterpart of such a fluctua- 
tion is an excess of positive spins on a transverse band of 
~ L sites in front of the approaching cluster, which has 
thus negligible probability when L — > oo. Starting from a 
localized cluster, the system spreads at constant speed be- 
cause the front shock is initially faster than the rear shock 
(see Fig. Hb, their velocity becoming equal only when the 
density in the band is uniformly equal to p h . Since the 
length of a high-density region is oc L, so does the time 
needed to reach the uniform profile. However, the mean 
reversal time t r {L) of a cluster does not increase propor- 
tionally to L in Id and clusters thus never fully spread in 
large enough systems: the total magnetization keeps flip- 
ping between m ~ ±m . 

We do not discuss the precise scaling of t r (L) BP since 
it turns out to be irrelevant for P(m). Indeed, while the 
time spent in a given cluster is proportional to t r (L), so is 
the time spent during a reversal. As shown in Fig. [3] a re- 
versal occurs because a fluctuation at the front of a cluster 
progressively flips all its sites. It ends when this fluctuation 
has gone through the whole cluster, whose length is pro- 
portional to t r (L). The whole process thus takes a times 
oc t r {L). In the large size limit, the ratio of the probability 
of finding the system in a cluster or in a reversal is thus con- 
stant, P(m) has a non-vanishing flat part between ±m 
(see Fig. Ell and (to) = 0: there is no spontaneous symme- 
try breaking in Id. A consequence is that the susceptibility 

X m = ^LK) ~ i 171 *) 2 = ^E 4 ( m ?) is extensive in 
the coexistence region, as can be seen in Fig. |4] 

Surprizingly, the finite-size scaling turns out to be more 
difficult in Id than in 2d. First, one cannot use |m, | in- 
stead of rrii when computing Xm, as is frequently done for 
the Ising model. The reversals are indeed important since 
they capture a finite part of the steady-state measure. Fur- 
thermore, since the magnetization varies continuously in 
the large L limit (which is easily reached in Id), the three 
peaks in P(m) can be very difficult to discriminate. When 
the width of the peaks becomes larger than their separation, 
no negative peak in G is observed. In figure|5]we show two 
extreme cases: without the RMFM to analyse the data, it 
would be very difficult to realize that one is looking at the 
same transition. 
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FIG. 5. Histograms and Binder cumulant of the total magneti- 
zation for po = 3, D = 1 (left) and p = 0.2, D = 10 (right). 
e = 0.9. L — 8000 for P(m). 

Another interesting difference between Id and 2d is that 
the m = to homogeneous profile is not stable in the ther- 
modynamic limit. Indeed, even though a fluctuation that 
would create a small negative cluster from a uniform pro- 
file to = to > is rare, its probability does not decay 
exponentially fast with L since one only has to flip a fi- 
nite number of sites. When the system size increases, the 
entropy of such a localized perturbation increases and the 
time it takes to exit the homogeneous state thus vanishes 
when L — > oo. In Id, there is thus only one transition line 
in the T, p Q ensemble, that separates homogeneous profiles 
from flipping clusters. 

Conclusion In this letter we have introduced a lattice 
model of self-propelled Ising spins, whose phenomenol- 
ogy is very similar to that of traditional models of aligning 
self-propelled particles. The simplicity of our model al- 
lows us to show that its flocking transition amounts to a 
liquid-gas phase transition with an infinite critical density. 
The total magnetization is proportional to the liquid frac- 
tion and thus varies continuously through this first-order 
phase transition, a rather counter-intuitive result. This sce- 
nario, confirmed by numerical simulations in 2d, is altered 
by the strength of fluctuations in Id, where no spontaneous 
symmetry breaking is observed. 

While the scenario presented here is consistent with the 
numerical results obtained on off-lattice systems, this is of 
course not a proof that these systems undergo the same type 
of phase-transition. Nevertheless, it already opens new 
perspective to study flocking models. For instance, giant- 
number fluctuations, which have been reported in many 
systems, are trivially present in the coexistence region of 
our model: since P{pi) is double-peaked, around p t and 
Ph, the variance of the number of particles in a box of 
finite size has to satisfy (N 2 ) - (N) 2 oc (A) 2 . Then, 
part of the difficulty of analysing the transition comes from 
the "choice" of the temperature-density ensemble where no 
discontinuous jump of the magnetization is seen in the ther- 
modynamic limit. The design of grand-canonical counter- 
parts of flocking models, in which the magnetization would 
jump discontinuously at a transition line, seems a promis- 
ing line of research, even though "changing ensemble" is 
not obvious out-of-equilibrium. 

The authors thank H. Chate, G. Gregoire, F. Peruani, F. 
van Wijland for useful discussions. 
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